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Abstract 

The cosmological models based on teleparallel gravity with nonzero torsion are considered. To 
investigate the evolution of this theory, we consider the phase-space analysis of the /(T) theory. It 
shows when the tension scalar can be written as an inverse function of x where x = pe/ {^m'^iH'^) 
and T = g{x), the system is an autonomous one. Furthermore, the We — phase analysis is 
given out. We perform the dynamical analysis for the models f{T) = /3Tln(T/To) and f{T) = 
'^"^p«(~-^/"^p«)" particularly. We find that the universe will settle into de-Sitter phase for both 
models. And we have examined the evolution behavior of the power law form in the Uep — <^ep 
plane. 
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I. INTRODUCTION 



Recent cosmic acceleration in our universe is suggested by a combination of different 
cosmic probes that primarily involves Supernova data . Besides the most well-known 
dark energy scenario (e.g., the scalar field quintessence model Q]) as the mechanism for late 
time cosmic acceleration, another kind of model based on infra-red modifications to general 
relativity (such as the f{R) theory 4-7|) was also proposed. The validity of general relativity 
on large astrophysical or cosmological scale has never been tested but only assumed, the 
recent cosmic acceleration in our universe might be nothing but the signal of breakdown of 
General Relativity at large scale js-ll]. 



The Riemann-Cartan geometry was proposed with the aim of unifying gravity and elec- 
tromagnetism 12| . In general, a large number of connections can be defined on a manifold. 



The assumptions of torsion-free and metric compatibility lead to the Einstein general rel- 
ativity with Levi-Civita connection. Meanwhile, the assumption of curvature-free leads to 
the teleparallel theory of gravity with Weitzenbock connection 
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16(. The /(T) theory 



is proved to be equivalent to general relativity but with different origins. The /(T) the- 
ory is considered as a gauge theory and general relativity is thought as a geometric theory. 
Recently, the /(T) gravity has be en p roposed to explain the present cosmic accelerating 



expansion without dark energy 



17H22]. Similar to f{R) gravity, based on a modification of 



the teleparallel equivalent of general relativity Lagrangian, the additional term related to 
the form of /(T) could be written as an effective energy. The torsion T in /(T) theory will 
be responsible for the observed acceleration of the universe. 

n 

The /(T) theory was firstly used ^_as the source of driving inflation 17]. Then it was 



applied to the late time acceleration 
was discussed in Refs. 
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2l|. Furthermore, the perturbation of /(T) theory 



2J], the reconstruction of the /(T) theory was presented in Refs. 



261] , and the observational constraint was analyzed in Refs. 27|, |28|] . In this letter, we 



made a dynamical analysis on two /(T) models by choosing variables different from those 



in Refs. 



19 



29|. 



This letter is organized as follows. In Sec. [TTl we introduce the /(T) theory, present 
the phase-space analysis and give out the analytical form of u'^. Two particular models are 
discussed in Sec. IIIII Finally, we draw our conclusion in Sec. HVl 
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II. THE /(T) GRAVITY 



In /(T) theory, the fundamental dynamical object is the vierbein field e^(x). In the 
teleparallelism, the orthonormal tetrad components e'X{x) are used, where the index A runs 
over 0, 1, 2, 3 in the tangent space at each point of the manifold, and fi is the coordinate 
index in the manifold and also runs over 0, 1, 2, 3. The curvature tensor and the covari- 
ant derivatives of e^(a;) with respect to the connection vanish globally, therefore e^(a;) are 



absolutely parallel vector fields, this theory is also called teleparallel gravity [l2j and the 
geometry is the Weitzenbock space-time characterized by torsion tensor alone. The vierbein 
field is related with the metric g^^ by 

gf,u = VABe^e^. (1) 
The torsion T^,,,, and contorsion K^'^ „ tensors are defined as 

= {d,et - d^e^) , (2) 
K^% = -i(T^-^-T'^%-T^g . (3) 

The action for /(T) theory is 

lT = ^jd'x\e\{T + f{T)), (4) 

where |e| = det {e^ = ^f—g-i the torsion scalar is 

^ = (5) 

and 

S^\^\ {k^'' , + 5; - 5; T"^) . (6) 
Varying the action with respect to the vierbein yields the equation of motion, 

+e-'d,{ee'^S^\){l + h) + c^S/^/tt^^T = 47rGe%T^ " (7) 

where T" is the energy-momentum tensor and the subscript T denotes the derivative with 
respect to the torsion scalar. 

For simplicity, we assume a flat Friedmann-Robertson- Walker metric, 

3 

ds"" = -dt" + a^{t)Y,{dx'f , (8) 



where a{t) is the scale factor. In the FRW space-time, Qf^^ = diag(— 1, a^, a^, a^) and therefore 
the tetrad components = (l,a,a, a) yield the torsion scalar T = — 6if^, where H = a/a 
is the Hubble parameter. 

The energy conservation equations for the radiation and pressureless matter are 

P; + 4p^ = 0, (9) 
p'^ + ?>pm = 0, (10) 

where and pm are the energy densities of radiation and pressureless matter respectively, 
and a prime means the derivative with respect to In a. From Eq. ([7]), the Friedmann equation 
for /(T) theory is gotten, 

= ^{Pm + Pr+ Pe). (H) 

where the effective (dark) energy density due to /(T) is 

Obviously, when / = 2T/t, i.e., /(T) oc T^/^, the effective energy density vanishes. The 
effective dark energy density is assumed to be conserved as well 

pl + 3(pe+Pe)=0, (13) 

where the effective pressure is Pe = —Pe — f^'^Afr + '^TfTT)T'/6, and the evolution of the 
torsion scalar T could be expressed as 

H = -1 3pm + 4:pr , . 

T 3mliH^2TfTT + fT + l' ^ ' 

A. The Phase-Space Analysis 

To make the phase-space analysis, we introduce the following dimensionless variable, 

3m2^i72' y-''"'- SmliH^' ' ~ ^mim' 

In terms of these dimensionless variables, the Friedmann equation becomes the constraint 
equation 

X + y + z=l. (16) 
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Combining Eqs. (|9|), (fTO|), (1T3|) and (jTi]), we get 

y' = -1/(3 + ^), 
^' = -^(4 + ^), 

and 

r _ {H^y _ 3y + 4z 

T ~ IP' ~ ~2TfTT+l + fT 

From the definition of x in Eq.f jTSj) . we know tliat once tfie function /(T) is specified, x can 
be expressed as a function of T , so if the inverse function exists, then T can be expressed 
function of x, 

T = g{x). (21) 

With the help of the above relation, Eqs. f|T7|) . f|T8|) and f|T9l) become an autonomous system. 
Because of the constraint ( fT6l) . there are only two independent variables, here we choose them 
as X and y. 



(17) 
(18) 
(19) 

(20) 



B. The Equation of State 



The elucidation of the physics of models becomes difficult when the EoS parameter is close 



to -1 



30|, since the evolving slowly effective energy density may be indistinguishable from 



ACDM scenario. By emphasizing the dynamics, the restricted regions of the trajectories 
are recovered in "position" and "velocity"- the value of the equation of state ratio u and its 
time variation u' 31\\ . 

The EoS parameter of the effective dark energy is, 

T' fx + '^T fxT 



= ^ = -1 + 

Pe ^ T3(-//r + 2/^) 



(22) 



When /(T) = constant, the model can be regarded as the cosmological constant theory 
which is attributed to the quantum zero-point energy of the particle physics vacuum. 

It is currently constrained from the distance measurements of SNIa as loq = —1.31'*'°'^^ 



and Uq = —lAt 



-0.90 
-0.81 



30|. In /(T) theories, the form of u'^ can be expressed as 



3{-f/T + 2fTrT 



3(-//r + 2/r) 3(-//T + 2/t)V T 



0.28 



(23) 



It is noticed that the parameters ui^ and ui'^ are functions of T, pm and p-y. But only 
the Priedmann equation could be used as a constraint equation. It is hard to draw the 
phase-space diagram since both the parameters ujf, and are functions of two variables. 
Fortunately, in the pivot periods, we could get the relation between — once the exact 
form of /(T) is given. In the radiation dominated phase, the energy of matter could be 
ignored; and once the matter is dominating, the energy density of radiation could be ignored. 
In the following , we will give two exact forms of f{T) as examples. 



III. TWO MODELS 



A. The Logarithmic Form 

Firstly, we consider a phenomenological model with a logarithmic form, 

/(r) = ^rin(^), (24) 
where ^ and Tq are constants. The effective energy density is 

and the dimensional variable x is 

Therefore, 

T = Toei'^+^ (27) 

and 



T' _ A-Ax-y 
Y ~ ~3/3 - 1 - xjP' 

The dynamical equations of the autonomous system become 



(28) 



«'--^(3- 3^:"^;^ )- (30) 

Setting x' — y' — 0, we get the critical point x — 1 and y = 0. The critical point corresponds 
to dark energy domination with — Q.^. — Q and We — —1. The eigenvalues are Ai = —3 
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and As = 4/37(3/32 - /3 - 1), so when (1 - ^)/6 < /3 < (1 + yT3)/6 and /3 0, the 
critical point is a stable fixed point. Unlike the standard model, there is no critical points 
for radiation domination z = 1 and matter domination ?/ = 1 for this model. In the view of 
dynamical analysis, this model is ruled out by the history of our universe. Therefore, it is 
not necessary to discuss its EoS parameter. 



B. The Power Law Form 



Now we consider the power law form, 

/(T) = am%{^r, (31) 

where a and n are dimensionless parameters. When n = 0, /(T) is a constant and the model 
becomes the ACDM model. By using the power-law form of /(T), we get the effective dark 
energy density, 

1 —T 

Pep = am%{n - 7^){—tT, (32) 
^ ''''pi 



and the effective pressure, 

Pep = -ampi{n--){—^) + ^(-^) T. (33) 



''pi " "V 

When n = 1/2, pep = Pep = 0. When n = 1, pep = = Samp^iY^, so the effective 

energy tracks the background matter. Substituting Eq. (152]) into the definition of the 
dimensionless variable x, we get 

x = a(2n-l)(^r-\ (34) 
^pi 

and 

/ a; \ ^/^""^^ 



So 

T' _ 4 - 4x - ?/ 
T 1 — nx ' 

and the dynamical equations are. 



(36) 



^ = l^x ^' 

y' = -y(3- \^^~^ ). (38) 

1 — 



Phases 


{x,y) 


Stability 


(Ai,A2) 


R 


(0,0) 


unstable 


(1, -8(n2 - 1)) 


M 


(0,1) 


unstable when n? < I 


(-1, 6(1 -n2)) 


A 


(1,0) 


stable 


(-3, -8) 



TABLE I: The properties of the critical points for the power law form with n ^ 1. 




FIG. 1: The phase diagram of attractor A for power law f(T) form with different initial conditions. 

Setting x' = y' = 0, we obtain the critical pints {x, y) = (0, 0), (x, y) = (0, 1) and 
(x, y) = (1, 0). The properties of the fixed points R, M, and A are summarized in Tabled 



These fixed points are similar to those found in Ref. 
be stable when > 1. 



19| . but we find the fixed point M can 



R The radiation dominated phase 

For the critical point R, x = y = and fi^ = 1. Because Ai = 1 > 0, it is an unstable 
fixed point, so the universe could exit from the radiation dominating phase. 

M The matter dominated phase 

For the critical point M, x = and y = 1. When > 1, it is a stable attractor, 
otherwise it is an unstable fixed point. So when < 1, the universe could exit from 
the matter domination. 



A The accelerating phase 
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-1.6 -1.4 -1.2 -1 

Wep 



-0.8 -0.6 -0.4 



FIG. 2: The trajectories of Wep — uj'^p in the phase space. From top to bottom, n takes the values 
of 0.3, 0.1, -0.1, -0.3, respectively. 

For the fixed point A, n 1, x = 1 and y = and it is a stable attractor independent 
of the parameters a and n. So Qm = = 0, We = —1 since T'/T = 0, and this phase 
recovers the de-Sitter phase. Now we have two attractors M and A when ra^ > 1. For 
a given initial condition, which attractor will the universe choose in the phase space? 
To answer this question, we go back to the dynamical Eqs. fl37|l and fl38|) . Note that 
there is a singularity in the system when x = 1/n. So if we start with x > 1/n, then 
the attractor is A; otherwise if we start with x < the attractor is M. When 
ra^ < 1, we have only one attractor A. We plot the phase diagram for the attractor A 
with n = 0.1 in Fig. [TJ 

For the special case n = 1, if a < 1, x = a. The constraint equation becomes y+z = 1 — a, 
so there is only one independent variable. And there are only two critical points R with 
= 1 — a and M with = 1 — a. 

Based on the power law form of f{T), the corresponding equation of state is obtained 

riT' 

Uep = -1 - - — . (39) 

In the radiation or the matter dominated phase, the evolution of the scale factor a is known, 
then it is easily to get the effective EoS parameter which is listed in Table [Tll Generally, 



Phases 


Scale factor 


T'/T 


Wei 




R 


a oc t^/^ 


-4 


-l + 4n/3 





M 


a oc 


-3 


-1 + n 






TABLE II: we summarize the critical points of the power law form. 



if the effective energy density part takes the dominated part, the corresponding effective 
EoS parameter Uei < which suggests n < 1. If the effective energy density part has 
the possibiUty to make our universe accelerate, the corresponding effective EoS parameter 
LOep < —1/3 which suggests n < 2/3. If the effective EoS parameter cuep < — 1, n < is 
suggested. 

After the universe transits to the matter dominated phase, the radiation part could be 
ignored, and the energy density of matter could be approximately written as 



Pn 



pi 



-T + f-2TfT). 



(40) 



Putting the above equations into Eg. (122!) . the equation of state is expressed as 

1 



n 



111 ~ 

^ep 



1 



riT 
3 Y' 



1 - an{2n - l)(-T/m2j"-i 
The derivative of the EoS parameter with respect to In a is changed to 



^ep — 



3(n- 1)^ 



n 



;i + Wep)( 



n 



I). 



CO. 



(41) 



ep 



After fixing the value of n, the diagram of Uep — oo'^p which is Figure [2] could be drawn out. 



There is a fixed point u, 



ep 



1, = which is corresponding to the accelerating phase A. 



IV. CONCLUSION 

In this letter, the phase-phase analysis of the /(T) theories is done and its stability of 
the critical points has been investigated. We can see that when the energy density can be 
written as an inverse function of x which could be expressed as a general form T = g{x), 
we can do the phase-analysis in the dynamical system. Furthermore, the value of the EoS 
parameter and its evolution are considered. Those analyses depend on the exact forms of 
/(T) heavily. 
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Specially, we made the dynamical analysis for the logarithmic form and the power form 
of f{T). For the logarithmic form /(T) = /3Tln(T/To), only one critical point which cor- 
responds to de-Sitter phase exists. Unlike the standard cosmology with Einstein gravity, 
there is no critical point for radiation and matter dominations. For the power law form 
/(T) = ampi{—T / rripi)"' , we find that two stable fixed points exist when > 1. One fixed 
point corresponds to matter domination and the other fixed point corresponds to de-Sitter 
phase. The reason for the existence of two stable fixed points is that the autonomous sys- 
tem has a singular point x = 1/n. If the system starts with x > 1/n, then the system 
will settle into de-Sitter phase A, otherwise the system will settle into matter domination 
M. When n = 1, de-Sitter phase is absent and the effective energy tracks the background 
matter. When ra^ < 1, only de-Sitter phase is the stable fixed point. And we have plotted 
the diagram of Uep — co'^p plane for the power law form. The fixed point Uep = —l,uj'^p = 
is corresponding to the accelerating phase A. 
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